Search for: All records

Abstract We describe internal, low‐frequency variability in a 21‐year simulation with a cloud‐resolving model. The model domain is the length of the equatorial Pacific and includes a slab ocean, which permits coherent cycles of sea surface temperature (SST), atmospheric convection, and the convectively coupled circulation. The warming phase of the cycle is associated with near‐uniform SST, less organized convection, and sparse low cloud cover, while the cooling phase exhibits strong SST gradients, highly organized convection, and enhanced low cloudiness. Both phases are quasi‐stable but, on long timescales, are ultimately susceptible to instabilities resulting in rapid phase transitions. The internal cycle is leveraged to understand the factors controlling the strength and structure of the tropical overturning circulation and the stratification of the tropical troposphere. The overturning circulation is strongly modulated by convective organization, with SST playing a lesser role. When convection is highly organized, the circulation is weaker and more bottom‐heavy. Alternatively, tropospheric stratification depends on both convective organization and SST, depending on the vertical level. SST‐driven variability dominates aloft while organization‐driven variability dominates at lower levels. A similar pattern is found in ERA5 reanalysis of the equatorial Pacific. The relationship between convective organization and stratification is explicated using a simple entraining plume model. The results highlight the importance of convective organization for tropical variability and lay a foundation for future work using coupled, idealized models that explicitly resolve convection. 

Abstract Cloud condensation and hydrometeor evaporation fractionate stable isotopes of water, enriching liquid with heavy isotopes; whereupon updrafts, downdrafts, and rain vertically redistribute water and its isotopes in the lower troposphere. These vertical water fluxes through the marine boundary layer affect low cloud climate feedback and, combined with isotope fractionation, are hypothesized to explain the depletion of tropical precipitation at higher precipitation rates known as the “amount effect.” Here, an efficient and numerically stable quasi‐analytical model simulates the evaporation of raindrops and enrichment of their isotope composition. It is applied to a drop size distribution and subcloud environment representative of Atlantic trade cumulus clouds. Idealized physics experiments artificially zero out selected processes to discern the separate effects on the isotope ratio of raindrops, of exchange with the environment, evaporation, and kinetic molecular diffusion. A parameterization of size‐dependent molecular and eddy diffusion is formulated that enriches raindrops much more strongly (+5‰ for deuterated water [HDO] and +3.5‰ for O) than equilibrium evaporation as they become smaller than 1 mm. The effect on evaporated vapor is also assessed. Rain evaporation enriches subcloud vapor by +12‰ per mm rain (for HDO), explaining observations of enriched vapor in cold pools sourced by evaporatively cooled downdrafts. Drops smaller than 0.5 mm evaporate completely before falling 700 m in typical subtropical marine boundary layer conditions. The early and complete evaporation of these smaller drops in the rain size distribution enriches the vapor produced by rain evaporation. 

Abstract High‐Resolution Multi‐scale Modeling Frameworks (HR)—global climate models that embed separate, convection‐resolving models with high enough resolution to resolve boundary layer eddies—have exciting potential for investigating low cloud feedback dynamics due to reduced parameterization and ability for multidecadal throughput on modern computing hardware. However low clouds in past HR have suffered a stubborn problem of over‐entrainment due to an uncontrolled source of mixing across the marine subtropical inversion manifesting as stratocumulus dim biases in present‐day climate, limiting their scientific utility. We report new results showing that this over‐entrainment can be partly offset by using hyperviscosity and cloud droplet sedimentation. Hyperviscosity damps small‐scale momentum fluctuations associated with the formulation of the momentum solver of the embedded large eddy simulation. By considering the sedimentation process adjacent to default one‐moment microphysics in HR, condensed phase particles can be removed from the entrainment zone, which further reduces entrainment efficiency. The result is an HR that can produce more low clouds with a higher liquid water path and a reduced stratocumulus dim bias. Associated improvements in the explicitly simulated sub‐cloud eddy spectrum are observed. We report these sensitivities in multi‐week tests and then explore their operational potential alongside microphysical retuning in decadal simulations at operational 1.5° exterior resolution. The result is a new HR having desired improvements in the baseline present‐day low cloud climatology, and a reduced global mean bias and root mean squared error of absorbed shortwave radiation. We suggest it should be promising for examining low cloud feedbacks with minimal approximation. 

Abstract The extension of a cloud‐resolving model, the System for Atmospheric Modeling (SAM), to global domains is described. The resulting global model, gSAM, is formulated on a latitude‐longitude grid. It uses an anelastic dynamical core with a single reference profile (as in SAM), but its governing equations differ somewhat from other anelastic models. For quasihydrostatic flows, they are isomorphic to the primitive equations (PE) in pressure coordinates but with the globally uniform reference pressure playing the role of actual pressure. As a result, gSAM can exactly maintain steady zonally symmetric baroclinic flows that have been specified in pressure coordinates, produces accurate simulations when initialized or nudged with global reanalyses, and has a natural energy conservation equation despite the drawbacks of using the anelastic system to model global scales. gSAM employs a novel treatment of topography using a type of immersed boundary method, the Quasi‐Solid Body Method, where the instantaneous flow velocity is forced to stagnate in grid cells inside a prescribed terrain. The results of several standard tests designed to evaluate the accuracy of global models with and without topography as well as results from real Earth simulations are presented.